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DIURNAL AND SEMIDIURNAL VARIATIONS IN EARTH ROTATION 
T. Herring 

Massachusetts Institute of Technology, Cambridge, MA 02139 


ABSTRACT 

During the last decade there has been an unprecedented improvement in both the accuracy and the temporal 
resolution of Earth rotation measurements. Determination of the position of the Earth's rotation axis both in 
inertial space and with respect to the crust with accuracies of about 0.3 milliarcseconds (mas) are now routine. In 
recent years, there has been and emphasis on the determination of short-period (daily and less) variations in Earth 
rotation. Two space based geodetic systems, very long baseline interferometry (VLBI) and the global positioning 
system (GPS) have proved to be very successful in this endeavor. Results for the tidally coherent part of the 
subdaily Earth rotation variations determined from the analysis of VLBI data are discussed. The magnitude of 
other subdaily variations are also considered, 

INTRODUCTION 

Variations in the rotation of the Earth are caused by the response of the Earth to external torques from the Sun, 
Moon and planets; by the response of the Earth to internal mass redistributions from surface and internal loads; 
and by exchanges of angular momentum between the solid Earth and its fluid elements (ocean, atmosphere and 
fluid inner core). (See /I / for a recent review.) 

In studies of the rotation of the Earth, two primary coordinates systems are used — one attached in specified 
manner to the surface of the Earth (the Terrestrial Reference System, TRS) and one attached to inertial space in 
which equations of motion of can be specified (the Celestial Reference System, CRS). Other intermediate 
coordinate systems are also useful and one attached to the rotation axis of the Earth is commonly used, although 
care must be taken when using this system to ensure the correct equations of motions are used and that the full 
transformation between the TRS and CRS is evaluated because invariably the rotation axis moves with respect to 
both of the systems. In all space geodetic analyses, it is necessary to evaluate the transformation between the TRS 
and CRS because the observations are made between objects whose coordinates change slowly or not at all in one 
system but not the other. For example, coordinates of Earth based receivers are slowly changing in the TRS 
whereas the coordinates of the transmitter (quasars or satellites) are either fixed or expressible in the CRS. To 
instantaneously transform between the TRS and CRS requires only three rotation angles be known assuming that 
the origins of the two systems are coincident. Conventionally, however, five angles have been used to express 
this transformation by introducing an intermediate coordinate system approximately aligned with the rotation axis 
of the Earth. The advantage of this latter approach has been that, at some level of accuracy, the angular 
arguments of the rotation matrices change slowly with time. (As we will see this is not the case when rotation 
variations of less than 1 milliarcsecond (mas) are studied). The conventional transformation between the TRS and 
CRS is given by: 


Rti = PNSVV (1) 

where R ti is the rotation matrix to instantaneously transform between the TRS and the CRS, and PNSW is the 
product of four rotation matrices. The \V matrix rotates the Z-axis of the TRS to the intermediate axis (called the 
Celestial Ephemeris Pole, CEP) using the X p and Y p components of polar motion, the S matrix rotates about the 
Z-axis of the coordinate system by an angle related to the rotation rate of the Earth and conventionally expressed 
in terms of Greenwich mean sidereal time (GMST) and the (small) difference between this rotation and the 
instantaneous direction of the Greenwich meridian as expressed by Universal Time; the PN matrix finally rotates 
the coordinate system to one which is fixed in inertial space using the nutation angles and precession. If the 
Earth's rotation was only affected by external torques whose frequencies when viewed from the TRS were all 
close to -1 cycles-per-sidereal-day (cpsd) and all other variations in the rotation occurred on time scales long 
compared to a day then all of the arguments of these rotation matrices would be changing slowly. However, this 
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is not the case when sub-milliarcsecond rotation variations are considered. Here we will use the conventional 
formulation as a starting point and express the more complicated rotational motion of the Earth as corrections to 
the angular arguments of the transformation in Equation (1). 

PARAMETERIZATION 

In conventional analyses, the argument of the S matrix, 0, (Greenwich sidereal time), is given by 

0 = GMSTohr + R UTl + Ay/cose (2) 

where GMST 0hr is Greenwich Mean Sidereal Time at 0 hours UTC, R is the ratio of sidereal time to universal 

time; UT1 is Universal time as measured by the rotation of the Earth, and the term Ay/cose is the equation of 
equinoxes /2, 3, 4/. Because only three rotation angles are needed to transform between the TRS and CRS, any 
rotational variations of the Earth can be represented by allowing continuous variations of the arguments Xp, Yp, 
and UT1. In practice, it is not possible to allow instantaneous, arbitrary variations because the space geodetic data 
depend on other parameters, such as atmospheric delay variations and clock variations, that also need to be 
estimated. However, by treating Earth rotation parameters as stochastic processes and enforcing some continuity 
of the values through the magnitude of the process noise variance, it is possible to estimate time varying Earth 
rotation parameters along with the parameters which need to be estimated (many of which are also stochastic 
processes /5, 6/. An example of the subdaily variations in Earth parameters obtained in such an analysis is shown 
in Figure 1. (The other results shown in Figure 1 will be discussed shortly). 

Analysis of a large number of VLBI experiments reveals that the subdaily variations in the Earth rotation 
parameters often appear dominated by diurnal and semidiurnal periodic signals. If this is the case, then the 
arguments of the W and S matrices can be written as 

UT\ = VT\ + wf, cos0 + «1, sin 0 + u c _ 2 cos2 G + u‘_ 2 sin 20 (3) 

X p = Y p - m c +1 cos 0 + m*, si n 0 - mf 2 cos 2 0 - ml 2 sin 2 0 - m\ 2 cos 2 0 + m' 2 sin 2 0 (4) 

Y p = Y p + m c +l cos0 + ml, sin 0 - ml 2 cos20 + ml 2 sin20 + m c +2 cos2G + m +2 sin 20 (5) 

where X~,Y~ and UTl are the slowly varying parts of the Earth rotation parameters and the coefficients ul,, ul,, 
uf 2 , and «1 2 are the cosine and sine components of the diurnal and semidiurnal variations in UTl; m‘, and m', 
are cosine and sine components of the prograde diurnal polar motion variations, m_ 2 and m_ 2 are cosine and sine 

components of the retrograde semidiurnal polar motion, and m\ 2 and m‘ +2 are the cosine and sine components of 
the prograde semidiurnal polar motion. The retrograde polar motion components have been retained as 
corrections to the nutation series in our analysis. Although, the coefficients in equations 3—5 are estimated with 
frequencies of exactly one and two cycles per sidereal day, the frequency range over which the admittance of the 
estimates is above 90% is very wide for 24 hour sessions of geodetic data and across the tidal bands (±0.1 cpsd 
about 1 and 2 cpsd) is greater than 98%. Specifically, a coherent signal with frequency/ will alias into these 

results with frequency /-f2 or f-2C2. depending upon whether f is in the diurnal or semidiurnal band. In Figure 2, 

we show estimates of u c _ 2 and u s _ 2 from the analysis 1085 VLBI experiments conducted between January 1984 
and June 1992. Although these estimates look largely like noise, a spectral analysis of the results (Figure 3) 
reveals that the results shown in Figure 2 are largely consistent with the noise in the estimates and are dominated 

by three spectral peaks which, to within the frequency resolution of the spectrum (3.3xl(H cpsd), correspond with 
the major semidiurnal tides. In Figures 4—8, we show the spectra of all of the diurnal and semidiurnal variations 
in Earth rotation parameters. In all cases the significant peaks correspond to major tidal lines and their relative 
amplitudes are approximately consistent with the tide generating forces, although in some cases there are some 
long period components which do not correspond to tidal lines. 

Such affects in Earth rotation were predicted and estimates made of the amplitudes in 1983 /7/ for the M2 tide, 
and by Brosche /8/ for three tidal constituents in each of the diurnal and semidiurnal bands. Given that these 
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variations appear to be tidally coherent, the parameterization of the diurnal and semidiurnal variations can be 
formulated such that only variations at the tidal frequencies are allowed. We have formulated the tidal arguments 
in the same fashion as those used for the nutation series with the explicit inclusion of a GST argument to account 

for the diurnal and semidiurnal variations. (For consistency with previous tidal work we use GST+Jt as the 
argument.) A general tide argument is then written as 

Pi = a A + b t £ +c i F + dp + e£l + f^O+K) ^ 

where a,. c,-, and /, are integer constants which multiply the conventional time dependent arguments of the 

motions of the Moon and Sun. The diurnal and semidiurnal variations in the rotation of the Earth can then be 
expressed as 


N 

8UT\ = ^u‘ cos <p i +u* sin <p i 
5X p = cos (p, + m\ sin (p k 

i = 1 

iV 

8Y p = ^ m‘ sin <p, + m\ cos <p, 

i=i 


(7) 

( 8 ) 

(9) 


where SUT1, SXp and 8Yp are the tidally coherent variations in UT1 and polar motion, the index i loops over N 
terms in both the diurnal and semidiurnal bands and the coefficients are quadrature components of each of the N 
tidal lines. (These coefficients should not be confused with the summed amplitudes in each of the tidal bands, 
equations 3-5). 

RESULTS 

To determine the amplitudes of the tidally coherent variations in the Earth’s rotation we have analyzed 1085 
VLBI experiments, each of approximately 24-hours durations, carried out between January, 1984 and June, 1992. 
On average, the experiments were separated by a few days and while this does mean that the full spectrum across 
the diurnal and semidiurnal bands can not be computed, it has little effect on determining the tidally coherent 
amplitudes because for these the frequencies are known. Within each of the tidal bands (including the prograde 
bands) we estimated the complex amplitudes at the eleven largest tidal lines except for the retrograde diurnal band 
where we estimated signals at the fourteen largest lines and at the frequencies corresponding to retrograde free- 
core nutation. The details of the type of nutational analysis are discussed in 19, 10,1 1/. In Tables 1 and 2 we give 
the values of these tidal contributions at the four major tidal lines in each of the tidal bands. We estimate that the 
standard deviation of these estimates arising from noise in the analysis of the VLBI data; from quasi-random high 
frequency variation in atmospheric angular momentum; and errors in the adopted models for the tidal 
displacements of the sites (solid-earth tides and horizontal and radial ocean loading) to be 1.4- 1.8 micro-time- 

seconds (pts) for the UT1 components and 20 micro-arc-seconds (pas) for the polar motion. These standard 
deviations correspond to site displacements of 0.6 mm. 

DISCUSSION 

The results presented in Table 1 are reasonably consistent with those presented in /1 2/ with the largest of the 
differences being attributed to (1) an increase in the duration of the VLBI data analyzed here and the removal of 
all VLBI experiments for which the rotational parameters could not be determined to better than one 
milliarcsecond; (2) the new analysis simultaneously estimates all of the coefficients accounting for the 
correlations, within a day session, between UT1 and polar motion coefficient estimates (the analysis in /1 2/ 
estimated the coefficients separately within each band from the time series within that band); and (3) our inclusion 
in this new analysis of corrections, in global sense, to the tidal displacement model used in the VLBI analysis. 
The largest of the corrections to the tidal model expressed as the amplitude of the appropriate (unnormalized) 
spherical harmonic were 1.5 mm radial and 0.3 mm horizontal, in the diurnal band (displacement at 45° latitude of 
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^ an ^ ^ mm rac ^ a ^ an ^ ^-2 mm horizontal in the semidiurnal band (displacement at 45° latitude 

of 0-8 “*0 °- 4 mm )- This type of correction to the tide model can accommodate global tidal model errors (such as 
the effects of mantle anelasticity), but can not account for errors in ocean loading models for example. This latter 
class of tidal model errors is difficult to account for because if all stations are allowed to have arbitrary tidal 
displacements, then the rotational component becomes undefined. (Any rotation variations could be represented 
by displacements of the sites.) One way of resolving this singularity is to adopt as the rotational component that 
which minimizes the tidal displacements at a selected set of sites. (The rotational component would then depend 
on the selected sites). Y 

For UT1, there are published values of the expected ocean tide contributions (Table 1) and our results are 
generally smaller than the ocean tide model values. Comparison between the observed rotation variations and 
ocean tide mcxlels will allow global constraints to be placed on ocean tide models and should provide information 
about some of the ad hoc parameters used in ocean tide modeling (see /8 /). Our results agree, with the standard 
deviations, of the results given in /13/, except at the K1 frequency where /1 3/ obtain a much larger estimate. 
I neir analysis used a much smaller data set and did not simultaneously estimate polar motion variations and tidal 
displacements. Sub-daily variations in Earth rotation parameters have also been observed using the GPS /JPL 
refs/, although at this time these analyses have only used a stochastic approach. 

In absolute terms the largest short period Earth rotation variations occur in the retrograde diurnal band because in 
this band the torques directly affect the rotation of the solid Earth. The response of the Earth in this band also 
strongly affects the rotational variations through the presence of the fluid core and the elastic and anelastic 
properties of the mantle (see /9, 10, 11/ for recent discussions). The only reason for tidally coherent torques 
outside of this band is because of the non-linear response of the oceans and atmosphere to applied torques at these 
high frequencies. 

While the subdaily variations, outside of the retrograde polar motion band, appear dominated by ocean tide 
effects, other sources of excitation also exist. Part of the M2 UT1 variation arises from libration in the Earth's 

rotation due to differences between the equatorial moments of inertia (1.9 pts /16 /), and 5 pas of the Ol prograde 
polar motion also arises from this same source /1 6/. There is also a contribution across all of these frequency 
bands from variations in atmospheric angular momentum. To quantitatively assess the magnitude of this 
contribution to the polar motion variations, we have integrated the polar motion excitations derived from 2 years 
of the twice-per-day estimates of AAM from National Meteorological Center (NMC) and spectrally analyzed the 
results. (The polar motion estimates were high-pass filtered with a Gaussian filter with a full- width-half 
maximum (FWHM) power width of about 60 days). The PSD function obtained from this analysis is shown in 
Figure 9 with spectra from the VLBI results superimposed. The excitation function included the wind 
contribution and assumed that the inverted barometer approximation was still valid at these high frequencies. 
(Assuming a rigid ocean increased the power in the high frequency region by approximately 50%. ) We may 
use this spectra to obtain estimates of the atmospheric contributions to high frequency polar motion. Integrating 
the power in the AAM spectrum for all frequencies corresponding to periods between 1 and 5 days, the AAM 

results would indicate that the standard of variations in polar motion with periods less than 5 days is 240 pas. 
This estimate includes the contributions from subdaily variations as well since their power will have been folded 
into the frequency range -1 to 1 cpsd. If the inverted barometer approximation is invalid at these high frequencies 

then the power in this frequency range would be 290 pas. (These estimates assume that NMC analysis has not 
smoothed the excitations excessively. The excitation spectrum show a flattening of spectrum for periods less than 
5 days suggesting that any smoothing of the results is probably small, and possibly indicating that much of the 
excitation may arise from noise in the NMC analysis.) The AAM spectrum can also be used to assess the 
atmospheric contribution to noise in the estimates of the diurnal and semidiurnal cosine and sine components. 
Ihe admittance bandwidth for these coefficients is approximately 0.8 cpsd (half power point), and the totai 

variance in a 0.8 cpsd bandwidth about 1 cpsd is 90 pas (108 pas with no inverted barometer correction), 
indicating that this much noise could enter into the estimates of the cosine and sine coefficients from a single 
VLBI experiment. We can not directly assess the contribution in the semidiurnal band, but it could not be larger 
than that tor the diurnal band because of the power folding discussed earlier. This is a non-negligible contribution 

to the best of the VLBI determinations which have standard deviations as low as 60 pas. It is also similar to the 
differences between daily estimates and the combined analysis results shown in Figure 1 . 
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CONCLUSIONS 

‘ h fJ° ng - hl | t0 7 of Y LB j experiments extremely accurate estimates of the tidally coherent variations of 
barth rotation in the diurnal and semidiurnal bands can be determined with the largest uncertainty in these results 
ansing trom the potential errors in the tidal displacements models used at the (small number) of VLBI sites. The 

VLBI results also show that remaining variation after removing the tidal signal is probably less than 300 pas (one 

standard deviation) for pole position and 20 pts for UT1. (Because of the wide admittance bandwidth of the 

STS S the dlU ™ al Se F“ si f nals ’ nearl y M of the P° wer of signals with frequencies greater than 1- 
L J e . pre ® cnt a ? the n ° ise level in the VLBI spectra.) However, in the period range between 1 and a few 
fmere^tinJ 1 !Ln!i con [ lt J u ? usl y sampled by VLBI because of the spacing of the experiments, there are likely to be 

denloSnt of ^S?i roc 10 ! he < ? ceanic . r , es P onse ‘o atmospheric pressure changes. In addition, the current 
p °yment of global GPS networks provides great promise for resolving more detail in the sub-daily spectral 

nf g TSr?c OU H h ‘w !L 0t clear , cur ntl y c th f the AAM estimates will be accurate enough to allow an interpretation 
of these results. With several years of data, GPS should also provide robust estimates of the tidally coherent 

v f r J at 1 K 7) s a * wed (except in the retrograde diurnal band) and will have the advantage over VLBI of far greater 
global distribution of stations. b 

wa i ^ up P° r !® d b y NASA Grant NAGW-0037, NOAA Grant NA90AA-D-AC481 , 
the Nor Grant EAR-89-05560, and by the Kerr-McGee Foundation. 
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TABIE L Observed and Theoretical Amplitudes of the Diurnal and Semidiurnal UT1 Variations. 


Fundamental Arguments 
It F D n o+n 

Tide 

Observed 
cos sin 

(pts) (gits) 

Brosch§ 
cos sin 

(pts) (gits) 

Diurnal UT1* 

1 0 2 0 2 -1 

01 

-17.3 

-16.1 

-14.2 

-36.3 

0 0 2 -2 2 -1 

PI 

-3.9 

-5.9 

-6.5 

-3.6 

0 0 0 0 0 -1 

K1 

6.5 

17.8 

15.0 

12.0 

1 0 2 0 2 -1 

Qi 

- 3.1 

-4.4 

Semidiurnal UTlt 

1 0 2 0 2-2 

M2 

-10.8 

14.3 

-20.0 

26.7 

0 0 2 -2 2 -2 

S2 

-0 1 

8.6 

-1.2 

16.2 

0 0 0 0 0 -2 

K2 

0.8 

3.7 

1 0 2 0 2 -2 

N2 

-1.6 

2.8 

-4.6 

5.9 


* The root mean square (RMS) scatter of the other 14 spectral components estimated is 1.5 (its. We estimate the 
standard deviations of the terms given to be 1 .6 gits 

t The RMS scatter of the other 14 spectral components estimated is 0.4 gits. We estimate the standard deviations 
of the terms given to be 1.4 gits. 

U ig n al s abo iif ? 39^ ^ S ca ' cu * at * ons account i n g for the effects of the fluid core which generally increases the 
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TABLE 2. Observed Amplitudes of the Diurnal and Semidiurnal Polar motion Variations. 


Fundamental Arguments 
t e F D FI 6+71 

Tide 

Retrograde^ 
cos sin 

(pas) (pas) 

Prograde§ 
cos sin 

(pas) (pas) 

Diurnal Polar Motion* 






0 0 

2 0 2 -1 

01 

-54 

-90 

-178 

89 

0 0 

2-2 2-1 

PI 

-38 

-7 

-49 

35 

0 0 

0 0 0 -1 

K1 



133 

-74 

1 0 

2 0 2 -1 

Ql 

14 

-24 

-33 

11 

Semidiurnal Polar Motion! 






0 0 

2 0 2 -2 

M2 

-10 

265 

1 

-58 

0 0 

2-2 2-2 

S2 

-67 

99 

2 

-12 

0 0 

0 0 0-2 

K2 

-26 

16 

39 

-5 

1 0 

2 0 2 -2 

N2 

-2 

11 

12 

-12 


! For the retrograde diurnal polar motion terms, the sin values correspond to in-phase nutation amplitude 
correction and the cosine terms to the out of phase nutation amplitude corrections. The apriori nutation series is 
ZMOA- 1990-2 discussed in /12/. 

§ For the prograde terms the signs of all of the fundamental arguments should be reversed. 

* The RMS scatter of the other 18 spectral components estimated in the retrograde diurnal band is 18 pas, and for 
the other 14 components in the prograde diurnal band 13 is pas. We estimate the standard deviations of the terms 
given to be 20 pas. 

t The RMS scatter of the other 14 spectral components estimated in the retrograde semidiurnal band is 8 pas, and 
in the prograde semidiurnal band is 8 is pas. We estimate the standard deviations of the terms given to be 20 pas. 
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Fig. 1. Subdaily variations in Earth rotation parameters obtained from the analysis of three VLBI experiments 
conducted the National Aeronautics and Space Administration (NASA) as part of the VLBI research and 
development program. In each of the panels, the dots with simple error bars are the results obtained by estimating 
(a) X p , (b) Yp, and (c) UT1 once-per-hour with the estimates treated as a random walk process with a process 
noise set such that adjacent values were constrained with standard deviation of 0.6 mas for X p , Yp, and 0.07 mts 
for UT1. The solid line shows the time variation predicted when the amplitudes of diurnal and semidiurnal 
variations are estimated instead of allowing stochastic variations. One standard deviation errors are shown 
occasionaly on these results. The dashed line shows the variations predicted by the tidally coherent model, given 
in Tables 1 and 2. 



21 

Oct, 1991 


22 


23 


24 




Sto (mas) Cos (mas) 


Herring 


10 


DIURNAL AND SEMIDIURNAL VARIATIONS IN EARTH ROTATION 


Fifi* 2. Estimates of cosine and sine components of the semidiurnal variations in UT1 from the analysis discussed 
in the text. One standard deviation error bars are shown. 
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Fig 3 Log of the Power-spectral density (PSD) function as a function of frequency for the semidiurnal UT1 
results shown in Figure 2. The thick (almost horizontal) line is the 99.5% confidence limit computed by 
propagating the variances in the data estimates through the weak smoother (Gaussian filter with full-width-at- 
half-maximum of 4 days) used to generate a uniformly spaced data set and through the fast Fourier transform used 
to generate the PSD. The correlations between the "smoothed" data were accounted for in computing the 
confidence interval. The spectral resolution is 1/3065 cpsd. For the UT1 spectra, we have expressed the UT1 
variations in angular arguments so that Figures 3 and 4 can be directly compared to Figures 5-8. 
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Fie 9 ComDOsite PSD function of polar motions variations for the frequency range ^ 

0.8 cpStendwidth is compressed into a 0.2 cpsd bandwidth m the VLBI spectra. 





